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We review recent theoretical progress aimed at understanding the formation and the early 
stages of evolution of giant planets, low-mass stars and brown dwarfs. Calculations coupling 
giant planet formation, within a modern version of the core accretion model that includes planet 
migration and disk evolution, and subsequent evolution yield consistent determinations of the 
planet structure and evolution. Uncertainties in the initial conditions, however, translate into 
large uncertainties in the luminosity at early stages. It is thus not possible to say whether young 
planets are faint or bright compared with low-mass young brown dwarfs. We review the effects 
of irradiation and evaporation on the evolution of short period planets and argue that substantial 
mass loss may have occurred for these objects. 

Concerning star formation, geometrical effects in protostar core collapse are examined by 
comparing ID and 3D calculations. Spherical collapse is shown to significantly overestimate 
the core inner density and temperature and thus to yield incorrect initial conditions for pre-main 
sequence or young brown dwarf evolution. Accretion is also shown to occur non-spherically 
over a very limited fraction of the protostar surface. Accretion affects the evolution of young 
brown dwarfs and yields more compact structures for a given mass and age, thus fainter 
luminosities, confirming previous studies for pre-main sequence stars. This can lead to severe 
misinterpretations of the mass and/or age of young accreting objects from their location in 
the HR diagram. Since accretion covers only a limited fraction of the protostar surface, we 
argue that newborn stars and brown dwarfs should appear rapidly over an extended area in 
the HR diagram, depending on their accretion history, rather than on a well defined birth line. 
Finally, we suggest that the distinction between planets and brown dwarfs be based on an 
observational diagnostic, reflecting the different formation mechanisms between these two 
distinct populations, rather than on an arbitrary, confusing definition. 



1. INTRODUCTION 

One of the fundamental questions of astrophysics re- 
mains the characterization of the formation of planets and 
stars. The mass ranges of the most massive planets and of 
the least massive brown dwarfs certainly overlap in the ^ 1- 
10 Mjup range; it is thus interesting to explore our under- 
standing of the planet and star formation mechanisms in a 
common review. 

The growing number of discovered extrasolar giant plan- 
ets, ranging now from neptune-mass to few Jupiter-mass 
objects, has questioned our understanding of planet forma- 
tion and evolution. The significant fraction of exoplanets 
in close orbit to their parent star, in particular, implies a 
revision of our standard scenario of planet formation. In- 



deed, these objects are located well within the so-called ice 
line and could not have formed in-situ. This strongly favors 
planet migration as a common process in planet formation. 
This issue is explored in §2 where we present consistent 
calculations between a revised version of the core accre- 
tion model, which does take planet migration into account, 
and subsequent evolution. In this section, we also review 
our current understanding of the effects of irradiation and 
evaporation on the evolution of short-period planets, hot- 
neptunes and hot-jupiters, and review present uncertainties 
in the determination of the evaporation rates. In §3, we 
briefly review our current understanding of protostellar core 
collapse and we show that non-spherical calculations are re- 
quired to get proper accretion histories, densities and ther- 
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mal profiles for the prestellar core. The effect of accretion 
on the early contracting phase of pre-main sequence stars 
and young brown dwarfs, and a review of observational de- 
terminations of accretion rates, are considered in §4. Fi- 
nally, through out this review, we have adopted as the defi- 
nition of planet an object formed by the three-step process 
described in §2.1, characterized by a central rocky /icy core 
built by accretion of planetesimals in a protostellar neb- 
ula. In contrast to genuine brown dwarfs, defined in this 
review as gaseous objects of similar composition as the par- 
ent cloud from which they formed by collapse. This issue 
is discussed in §5 and observational diagnostics to differ- 
entiate brown dwarfs from planets, based on their different 
formation mechanisms, are suggested. Section 6 is devoted 
to the conclusion. 

2. GASEOUS PLANETS: BIRTH AND EVOLUTION 



2.1. Planet formation 

The conventional planet formation model is the core ac- 
cretion model as developed by Pollack et al. (1996, here- 
after P96). One of the major difficulties faced by this model 
is the long timescale necessary to form a gaseous planet 
like Jupiter, a timescale significantly larger than typical 
disk lifetimes, < 10 Myr. Reasonable timescales can be 
achieved only at the expense of arbitrary assumptions, Uke 
e.g. nebula mean opacities reduced to 2% of the ISM value 
in some temperature range or sohd surface density signif- 
icantly larger than the minimum mass solar nebula value 
(Hubickyj et al, 2005). This leaves the standard core ac- 
cretion model in an uncomfortable situation. This model 
has been extended recendy by Alibert et al. (2004, 2005, 
hereafter A05) by including the effects of migration and 
disk evolution during the planet formation process. The oc- 
curence of migration during planet formation is supported 
by the discovery of numerous extrasolar giant planets at 
very short distance to their parent stars, well within the so- 
called ice Une, about 5 AU for the solar nebula conditions. 
Below this limit, above ice melting temperature, the insuffi- 
cient surface density of solids that will form eventually the 
planet core, and the lack of a large reservoir of gas prevent 
in-situ formation of large gaseous planets. 

Moreover, inward migration of the planet should arise 
from angular momentum transfer due to gravitational inter- 
actions between the gaseous disc and the growing planet 
(Lin and Papaloizou, 1986; Ward, 1997; Tanaka et al, 
2002). Taking into account the migration of a growing 
planet solves the long lasting timescale problem of the core 
accretion scenario. Indeed, when migration is included, the 
planet feeding zone never becomes depleted in planetesi- 
mals. As a result, the so-called phase 1 (see P96), domi- 
nated by accretion of sohd material, is lengthened whereas 



phase 2, dominated by gas accretion, is shortened appre- 
ciably. During the last so-called phase 3, runaway gas ac- 
cretion occurs and the predominantly H/He envelope is at- 
tracted onto the core. Phase 3 is very short compared to 
phases 1 and 2, and phase 2 essentially determines the for- 
mation timescale of the planet. The planet can thus form 
now on a timescale consistent with disk lifetimes, i.e. a few 
Myr for a Jupiter (see A05). 

In the models of Bodenheimer et al. (2000a) and Huby- 
ckij et al. (2005), which are based on the P96 formal- 
ism, the calculations proceed in 3 steps: (i) the planet is 
bounded by its Roche lobe {Rp = Rl) (or more precisely 
by Min(i?i, i?acc) where Race = GM/cl is the accretion 
radius and Cs the local sound velocity in the disk) so that 
the temperature and pressure at the planet surface are the 
ones of the surrounding nebula. Note that in P96 calcu- 
lations, opacity of the nebula is a key ingredient; (ii) the 
planet external radius is the one obtained when the max- 
imum gas accretion rate is reached. In P96, this value is 
fixed to 1 X 10~^ yr~^. At this stage, the external 
conditions have changed {Rp < Rl). Matter falls in free 
fall from the Roche lobe to the planet radius, producing a 
shock luminosity; (iii) once the planet reaches its prede- 
fined final mass, the accretion rate is set to and the bound- 
ary conditions become the ones of a cooling isolated object, 
L = Atto-R'^T^^^ and n^Pph = where kh denotes the 
mean Rosseland opacity. The planet surface radius is es- 
sentially fixed by the accretion shock conditions (see e.g. 
Fig. Id of Hubickyj et al., 2005). This value, however, re- 
mains highly uncertain, as its correct determination would 
imply a proper treatment of the radiative shock. In A05, 
phase (i) is similar to step (i) described above, except that 
the planet migration from an initial arbitrary location and 
the disk evolution are taken into account, so that the ther- 
modynamic conditions of the surrounding nebula, as well 
as the distance to the star, and thus the planet Roche lobe 
radius, change with time. The planet's final mass is set by 
the accretion rate limit, and is thus not defined a priori. Note 
that, because of the disk evolution and/or the creation of a 
gap around the planet, the accretion rate limit is 1 to 2 or- 
ders of magnitude smaller than the one in P96 at the end 
of phase (i) and reaches essentially with time, a fact sup- 
ported by 3D hydrodynamical simulations (D 'Angela et al. , 
2003; Kley and Dirksen, 2005). Eventually the planet opens 
a gap when its Hill radius becomes equal to the disk density 
scale heigth and migration stops or declines until the disk 
is dissipated (see A05 for details). The planet radius can- 
not be defined precisely in this model as it results from the 
competing effects of gas accretion and planet contraction 
with changing boundary conditions as the planet migrates 
inward and the disk evolves. In any event, the final stages of 
accretion are Ukely to occur within streams (see e.g. Lubow 
et al, 1999), i.e. non-spherically and, as mentioned above, 
the planet final radius remains highly uncertain, at least in 
any ID calculation. 

The migration rate, in particular type I migration for low- 
mass planet seeds, remains an ill-defined parameter in these 
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calculations. The observed frequency of extrasolar plan- 
ets implies a rate significantly smaller than estimates done 
for laminar disks {Tanaka et ai, 2002). Numerical mod- 
elling of turbulent disks yields significantly reduced migra- 
tion rates (Nelson and Papaloizou, 2004, see also D' Angela 
et at., 2003). It has been suggested recently that stochas- 
tic migration, i.e. protoplanets following a random walk 
through the disk due to gravitational interaction with turbu- 
lent density fluctuations in the disk, may provide a means 
of preventing at least some planetary cores from migrating 
into the central star due to type I migration (Nelson, 2005). 
Based on these arguments, and for lack of better determi- 
nations, AOS divide the aforementioned rate of Tanaka et 
al. by a factor 10 to 100. As noted by these authors, nu- 
merical tests show that, provided the rate is small enough 
to preserve planet survival, its exact value affects the extent 
of migration but not the formation timescale, nor the planet 
final structure and internal composition. 

2.2. Planet evolution 

2.2.7. Non irradiated planets 

We first examine the evolution of young planets far 
enough from their parent star for irradiation effects to be ne- 
glected. In order for the evolution to be consistent with the 
formation model, the planet structure includes now a cen- 
tral core surrounded by an envelope enriched in heavy ele- 
ments. These conditions are given by the formation model 
described in §2. 1, performed for different initial parameters 
(initial orbital distance, dust-to-gas ratio in the disc, photo- 
evaporation rate, disc initial surface mass). The planets are 
found to form with essentially the same core mass (Mcore — 
6 M^) independent of the planet final mass, whereas the 
heavy element mass fraction in the envelope deposited by 
the accreted planetesimals is found to increase substan- 
tially with decreasing total mass (Barajfe et al, 2006). The 
hydrogen-helium equation of state (EOS) is the Saumon, 
Chabrier and VanHorn EOS (Saumon et al, 1995) whereas 
the thermodynamic properties of the heavy material rele- 
vant to the planet structure (ice, dunite(= Mg2Si04), iron) 
are calculated with the ANEOS EOS (Thompson and Lau- 
son, 1972). In the present calculations, we assume that the 
core is made of dunite, as representative of rock, yielding 
typical mean densities in the core ~ 6-7 g cm""^. Compar- 
ative calculations with water ice cores, corresponding to a 
lower mean density ^ 3 g cm^"^, change only slightly the 
mass-radius relationship for planets of identical core and to- 
tal mass. As mentioned above, the specific heat of the core 
is calculated with the ANEOS EOS so that the core con- 
tributes to the planet thermal evolution. Fig. 1 displays the 
evolution of the radius and luminosity for 1 and 4 jupiter- 
mass planets, respectively. The solid and long-dash lines 
correspond to different initial radii for the new born planet, 
namely 3 and 1.3 Rj for the 1 Mjup planet and 4 and 1.3 



i?j for the 4 Mjup planet, respectively. The 1 .3 i?j case is 
similar to the calculations of Fortney et al. (2005), based 
on the aforementioned formation model of Hubickyj et al. 
(2005). Note that these values are comfortably smaller than 
the Roche lobe limits at 5.2 AU from a Sun (~ 530 i?j and 
~ 830 i?j for a 1 Mjup and a 4 M jup planet, respectively 
(Eggleton, 1983)). The t = age for the planet evolution 
corresponds to the end of its formation process, just after 
the runaway gas accretion (phase 3) has terminated. This 
planet formation timescale, namely ~2-3 Myr, should thus 
be added to the ages displayed in Fig. 1 for the planet evo- 
lution. As seen in the figure, the difference between these 
initial conditions, namely a factor ~2-3 in radius, affects 
the evolution of the planet for 10^ to 10^ yr, depending 
on its mass. This reflects the significantly different thermal 
timescales at the begining of the evolution (t = 0) for the 
different initial radii, namely Ikh = GA'P/RL = 3 x IQ-'' 
and ~ 5 X 10'' yr, respectively, for 1 Mjup. The smaller the 
initial radius the larger the consequences. Unfortunately, 
as mentioned above, uncertainties in the models of planet 
formation prevent an accurate determination of the initial 
radius of the new bom planet. Changing the maximum ac- 
cretion rate or the opacity in P96, for example, or resolving 
the radiation transfer in the accretion shock, will very likely 
affect the planet radius within a large factor. Therefore, at 
least within the present uncertainties of the planet forma- 
tion models, young gaseous planets with cores and heavy 
elements in their envelopes can easily be 10 times brighter 
than suggested by the calculations of Fortney et al. (2005) 
and thus are not necessarily "faint" in the sense that they 
can be as bright as pure gaseous, solar composition H/He 
objects of the same mass, i.e. low-mass brown dwarfs. In 
the same vein, the initial gravity of the planet can not be 
determined precisely and can certainly vary within at least 
an order of magnitude between log g ^ 2 and log g ^ 2) for 
a jupiter-mass. Detections of young exoplanet luminosities 
with reasonable age determinations, i.e. within < 10 Myr 
uncertainty, for instance in young clusters, would provide 
crucial information to help narrowing these uncertainties. 

2.2.2. Effect of irradiation 

We now examine the effects of irradiation on the evo- 
lution of close-in exoplanets, the so-called "hot-jupiters" 
and "hot-neptunes" objects. Inclusion of the effect of ir- 
radiation of the parent star on the structure and evolution 
of short period exoplanets has been considered by several 
authors. Only a few of these calculations, however, are 
based on consistent boundary conditions between the inter- 
nal structure and the irradiated atmosphere profiles. Such 
a proper boundary condition, implying consistent opaci- 
ties in the atmosphere and interior structure calculations, is 
determinant for correct evolutionary calculations of irradi- 
ated planets because of the growing external radiative zone 
which pushes the internal adiabat to deeper levels (Guillot et 
ah, 1996; Seager and Sasselov, 1998; Barman et al., 2001, 
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Fig. 1. — Evolution of the radius and the luminosity for a 1 Mjup (left) and a 4 Afjup (right) planet with a 6 solid core and 
A'/z,cnv/A^cnv = 10%, for two different initial radii (solid vs dash lines, see text). The dot-dash lines portray the cooling of coreless, pure 
gaseous brown dwarfs of solar composition with similar initial radii as for the solid lines; the differences reflect the influence of the 
presence of a central core on the evolution. 



2005). The out-going flux at the surface of the planet now 
includes the contribution from the incoming stellar flux Ti,\ 

R 

^out — CrT^ff + Tine = CrTc'ff + f{—)'^J^i, 

a 

= aT^g + {1 - A)J^,nc + AJ^.nc. (D 

In Eq.l, crT^ff denotes the intrinsic internal flux of the 
planet, A the Bond albedo and the last term on the r.h.s. 
of the equation is the reflected part of the spectrum. The 
factor / is a geometrical factor characterizing the stellar 
flux redistribution over the planet surface (/= 1 implies the 
flux is redistributed over vr steradians, /=l/2 that it is redis- 
tributed over the day-side only, as intuitively expected for 
tidally locked planets, and /=l/4 over the entire planet sur- 
face). Burkert et al. (2005) have performed hydrodynamic 
calculations related to the heating of the night side of syn- 
chronously locked planets. With reasonable assumptions 
for the opacity in the atmosphere, these authors find that the 
temperature difference between the day side and the night 
side could be in the ^ 200-300 K range, not enough to make 
an appreciable difference in the radius. Previous estimates 
{Showman and Guillot, 2002; Curtis and Showman, 2005; 
Iro, Bezard and Guillot, 2005), however, predict day/night 
temperature differences about twice this value, and this is- 
sue needs to be further explored. From Eq. 1, the evolution 
of the irradiated planet now reads: 

L=- T— + ^TTRlaT^g + Lreflected, (2) 



where = ^T^nc = (1 - A)f{^)^T^ denotes 
the planet equilibrium temperature, i.e. the temperature it 
would reach after exhaustion of all its internal heat content 
and contraction work (Tcff ^ 0). 

As shown in Chabrier et al. (2004) and Baraffe et al. 
(2005), consistent calculations between the iiTadiated at- 
mospheric structure and the internal structure, which fixes 
the boundary condition for the planet photospheric radius, 
reproduce the radii of all observed transit planets so far, 
without additional sources of internal heating, except for 
HD209458b, which remains a puzzle (see Fig. 1 of Baraffe 
et al, 2005). These calculations were based on planet in- 
terior models composed entirely of hydrogen and helium 
and do not include either a central core or heavy element 
enrichment in the envelope. The effect of a central rocky 
core on irradiated planet evolution has been examined by 
Bodenheimer et al. (2003) but with simplified (Eddington) 
boundary conditions between the atmosphere and the in- 
terior These authors found that for planets more massive 
than about 1 Afj^p the decrease in radius induced by the 
presence of a core is about 5%, in agreement with previous 
estimates for non-irradiated planets (Saumon et al, 1996). 
The effect, however, will be larger for less massive planets, 
including the recently discovered hot-Neptunes. This issue 
has been addressed recently by Baraffe et al. (2006), with 
proper, frequency-dependent atmosphere models. These 
authors find that, for a Saturn-mass planet (^ 100 M®), 
the difference in radius between a pure H/He planet and a 
planet with a 6 core and a mass fraction of heavy el- 
ement in the envelope Z=Mz_onv/A^nv=10%, as predicted 
by the formation model, is Rz /Rhhc — 0.92, i.e. a ~ 9% 
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effect, possibly within present Umits of detection. 

A point of concern in the present calculations is that the 
boundary condition between the irradiated atmospheric pro- 
file and the interior profile is based on atmosphere models 
of solar composition. Most of the transiting planets, how- 
ever, orbit stars that are enriched in metals and the planet 
atmosphere is supposed to have the same enrichment. Cal- 
culations including such an enrichment are under work (see 
§5). The effect, however, is Ukely to be small for two rea- 
sons. First of all, the enrichment of the parent stars remain 
modest, with a mean value [M/H\ «0.2-0.3 {Santos et ah, 
2004). Second of all, irradiated atmospheric profiles display 
an extensive radiative zone (see above) so that gravitational 
settling may occur even though, admitedly, various mixing 
mechanisms (e.g. decay of gravitational waves, convective 
overshooting, winds) could keep gaseous heavy elements 
suspended in radiative regions. Planets at large enough 
orbital distances for the effect of irradiation on the atmo- 
spheric thermal profile to be negligible, however, should 
display significant heavy element enrichment in their atmo- 
sphere, as observed for the giant planets of the solar system. 

2.2.3. Evaporation 

The question of the long-term stability of gaseous close- 
in extrasolar giant planets has been raised since the discov- 
ery of 5 1 Peg b. In the framework of Jeans approximation, 
the evaporation rate $ (hydogen atoms cm~^ s~^) is given 
by Chamberlain and Hunten (1987): 

$ = ^\[^^eM-X){l + X), (3) 

where riexo and Texo are the number density and the temper- 
ature at the exobase (the level at which the mean free path of 
hydrogen atoms equals the scale height) and X = / Vq 
is the escape parameter, Voo = {2GMp/ Rp)^/'^ the planet 
escape velocity and vq = (2fcT/m)^/^ the mean thermal 
velocity at Tcxo- The first estimates of the evaporation rate 
of hot Jupiters (e.g. Guillot et al, 1996) were obtained by 
using the equilibrium temperature Teq instead of the un- 
known value of Tcxo- For a typical 5 1 Peg-b-like hot Jupiter, 
(IMjup, Teq ~1300 K), the escape parameter X = vl^/vl 
is then found to be larger than 150 whereas escape rates 
become significant for values below 20. On this basis, hot 
Jupiters were claimed to be stable over the lifetime of their 
star. However, Teq is not the relevant temperature for ther- 
mal escape, which occurs in the exosphere, where heating is 
due to XUV irradiation. With simple assumptions, several 
authors estimated that the exospheric temperature could be 
of the order of 10,000 K (X < 20) and thus attempted the 
observation of the escaping H (Moutou et al., 2001). Lam- 
mer et al. (2003, L03) showed that the conditions allow- 
ing the use of Jeans approximation (hydrostatic equilibrium 
and neghgible cooling by the escape itself) are not met in 
hot Jupiters, because of the considerable heating by stellar 



XUV. The application of Jeans escape yield unrealistically 
high exospheric temperatures {X < 1) in contradiction 
with the required hydrostatic hypothesis. They concluded 
that hot Jupiters should experience hydrodynamic escape, 
without a defined exobase, where the upper atmosphere is 
continuously flowing to space and maintained at low tem- 
perature (^ 10, 000 K) by its expansion. In this blow-off 
model, the escape rate of the main atmospheric component, 
H, is only limited by the stellar XUV energy absorbed by 
the planet and is given by: 

M = 3(^^ye^,/(Gp), (4) 

where p is the mean planetary density and !Fi, is the stel- 
lar flux, averaged over the whole planet surface, including 
both the contribution in the 1-1000 A wavelength interval 
and the 1215 A Lyman-a line. Rxuv is the altitude of the 
(infinitely thin) layer where all the incoming XUV energy is 
absorbed while Rp is the radius observed in the visible dur- 
ing a transit. Here, e would represent the heating efficiency, 
or the fraction of the incoming XUV flux that is effectively 
used for the escape. L03 applied a hydrodynamic model 
(Watson et al, 1981) and estimated Rjcc/y/Rp ~ 3 for 
orbital distances closer than 0.1 AU. By assuming e = 1 
(or, in other words, that escape and expansion are the only 
cooling processes) they inferred the physical upper Umit 
for the XUV-induced thermal escape rate to be 10^^ g/s for 
HD209458 b at present time. Considering the evolution 
of XUV emission of main sequence G stars {Ribas et al., 
2005) and the significantly lower density of young gaseous 
planets implies rates 10 to 100 times higher in the early his- 
tory of the hot Jupiters. Using these simple arguments, L03 
suggested that hot Jupiters could have been initially much 
more massive although more detailed models are needed to 
better estimate the effective hydrodynamic escape rate. 

Independently of this theoretical approach, Vidal-Madjar 
et al. (2003, VM03) measured the absorption in the Lyman- 
a line of HD 209458, using STIS onboard HST, during the 
transit of its planet. The decrease of luminosity they found 
is equivalent to the transit of a RLya = 3 Rp opaque disk. 
Although this observation seems to be consistent with L03, 
a larger but optically thin hydrogen cloud can also account 
for the observation. In fact, by noticing that the Roche 
lobe radius of the planet was 3-4 Rp, VM03 concludes that 
part of the observed hydrogen must consist in an escap- 
ing cometary-like tail. They estimated that the absorption 
implies an escape rate not lower than 10^" g/s. 

The truncation of the expanded atmosphere by the Roche 
lobe, which was not considered by L03, has obviously to 
be taken into account in the mass loss process. Lecave- 
lier et al. (2004) proposed a geometrical blow-off model in 
which a hot exobase (~ 10, 000 K), defined according to 
Jeans approximation, reaches the Roche lobe radius. This 
yields enhanced loss rates compared to a classical Jeans cal- 
culation that would not take into account the gravity field 
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and the tidal distorsion of the atmosphere. Jaritz et al. 
(2005) argued that, although geometrical blow-off should 
occur for some of the known hot Jupiters, HD 209458b ex- 
pands hydrodynamically up to 3 Rp without reaching the 
LI Lagrange point at which the Roche lobe overflow oc- 
curs. If confirmed, the debated observation of O and C in 
the expanded atmosphere of HD 209548 b {Vidal-Madjar 
et al., 2004) would favor the hydrodynamic regime, which 
is required to drag heavy species up to the escaping lay- 
ers. However, the STIS instrument is no longer operational 
and similar observations will have to wait new EUV space 
observatories. Another indirect confirmation of the hydro- 
dynamical regime is the absence of an Hq, signature beyond 
Rp {Winn et al, 2004). This can be explained by the low 
temperature (< 5000 K) expected in the hydrodynamically 
expanding atmosphere. Yelle (2004) published a detailed 
model of the photochemistry, radiative budget and phys- 
ical structure of the expanding upper atmosphere of hot 
Jupiters and derived a loss rate of 10^ g/s, about a factor 
100 lower than the value inferred by VM03 from the obser- 
vation. Recendy, Tian et al. (2005) pubhshed an improved, 
multi-layer hydrodynamical model (compared to Watson), 
in which the energy deposition depth and the radiative cool- 
ing are taken into account. Rates of the order of 5 x 10^° g/s 
are found, although they also depend on an arbitrary heating 
efficiency e. It is important to note that the composition of 
the expanding atmosphere in heavy elements can dramati- 
cally affect its behavior, mainly by modifying the radiative 
transfer (absorption and cooling). 

Non-thermal escape is much more difficult to estimate 
as it depends on the unknown magnetic field of the planet 
and stellar wind. Thermal escape is usually considered as 
the dominant mass loss process {Griefimeier et al., 2004), 
but considering the complexity of the magnetic coupling 
between the star and the planet at orbital distances closer 
than 0.045 AU, unexpected non-thermal processes may still 
dominate the evaporation of some short-period exoplanets. 

VM03 and L03 both suggested that the evaporation 
could lead to the loss of a significant fraction of the ini- 
tial planetary mass and even to the evaporation of the whole 
planet, possibly leaving behind a dense core. In order to in- 
vestigate the possible effects on the mass-radius evolution 
of close-in exoplanets, Baraffe et al. (2004, 2005) included 
the maximum XUV-limited loss from L03 in the simulated 
evolution of a coreless gaseous giant planet, taking also into 
account the time dependency of the stellar XUV luminos- 
ity, calibrated on observations (Ribas et al, 2005). These 
studies showed that, even at the maximum loss rate, evap- 
oration affects the long-term evolution of the radius only 
below an initial critical mass. For initial masses below 
this critical mass, the planet eventually vanishes in a very 
short but dramatic runaway expansion. This critical mass 
depends of course on the escape rate considered and drops 
to values much below 1 Mjup when using lower rates hke 
the ones predicted by Yelle, Tian et al., and Lecavelier et al. 
(Baraffe et al, 2006). One interesting result of the Baraffe 
et al. work needing further attention is that evaporation 



does not seem to explain the surprisingly large visible ra- 
dius (Rp) of HD 209458b, except if this planet is presently 
seen in its last and brief agony, which seems extremely 
unlikely. The explanation for the large observed radius of 
HD 209458b thus remains an open question. 

One may wonder whether this runaway evaporation 
phase can be studied with hydrostatic atmosphere mod- 
els and quasi-static evolution models. Atmospheric hydro- 
static equilibrium is valid for values of the escape parameter 
X > 30. For a hot Jupiter at 0.045 AU, values of X below 
30 are found in the thermosphere, where the temperature 
increases above 7000 K, at i? > 1.1 Rp (see for instance 
Yelle, 2004). Such levels, with number densities n < 10® 
cm^*^, lie well above the levels where the boundary condi- 
tion applies, i.e. near the photosphere with gas pressures 
P ~ 10~^-10 bars. The quasi-static evolution assumption 
is justified by the fact that, even though the characteristic 
timescale of evaporation, M/M, can become comparable 
to or even shorter than the Kelvin-Hehnholtz timescale, 
^KH ~ IGrr? I (RL), it remains much larger than any hy- 
drodynamical timescale. The present runaway phase, in- 
deed, refers to a thermal runaway, hke e.g. thermal pulses 
in AGB stars, characterized by a thermal timescale. Quasi- 
static evolution thus remains appropriate to study this mass 
loss process, at least until truly hydrodynamic processes 
affect the planet photosphere. 

More recently, Baraffe et al. (2006) examined the pos- 
sibility for lower mass hot-neptune planets (1 Mnbp = 
18 M0 ^ 0.06Afj„p) to be formed originally as larger 
gaseous giants which experienced significant mass loss dur- 
ing their evolution. Depending on the value of the evapo- 
ration rates, these authors showed that presently observed 
(few gigayear old) neptune-mass irradiated planets may 
originate from objects of over a hundred earth masses if 
the evaporation rate reaches the maximum L03 value. For 
~ 10-20 times lower rates, as suggested e.g. by the hy- 
drodynamical calculations of Tian et al. (2005), the hot- 
Neptunes would originate from objects of ^ 50 M®, mean- 
ing that the planet has lost more than 2/3 of its original 
mass. For rates a factor 100 smaller than L03, the effect of 
evaporation is found to become more modest but a planet 
could still loose about 1/4 of its original mass due to stellar 
induced evaporation. These calculations, even though ham- 
pered by the large uncertainty in the evaporation rates, show 
that low-mass irradiated planets which lie below the afore 
mentioned critical initial mass, may have originally formed 
as objects with larger gaseous envelopes. This provides 
an alternative path to their formation besides other scenar- 
ios such as the core-colUsion model {Brunini and Cionco, 
2005). 

3. GRAVITATIONAL COLLAPSE OF PRESTEL- 
LAR CORES 

After having examined the status of planet formation and 
evolution, we now turn to the formation and the early stages 
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of evolution of stars and brown dwarfs. In this section, we 
first review our current knowledge of the gravitational col- 
lapse of a protostar. We then will focus on the importance 
of non-spherical effects in the collapse. 

3.1. One dimensional models 

Numerous authors have extensively considered the ID 
collapse of a spherical cloud. One of the most difficult as- 
pects of the problem is the treatment of the cooling of the 
gas due to colhsional excitation of gas molecules, partic- 
ularly during the late phase of the collapse when the gas 
becomes optically thick. Radiative transfer calculations 
coupled to hydrodynamics are then required. However, as 
noted originally by Hayashi and Nakano (1965) and con- 
firmed by various calculations {Larson, 1969; Masunaga 
and Inutsuka 1998; Lesaffre et ai, 2005) the gas remains 
nearly isothermal for densities up to 10^-10^ cm~^, mak- 
ing the isothermal assumption a fair and attractive simpUfi- 
cation. 

3.1.1. The isothermal phase 

The isothermal phase has been extensively investigated 
both numerically and analytically. In particular, a family 
of self-similar solutions of the gravitational contraction has 
been studied in detail by Pension (1969), Larson (1969), 
Hunter (1977), Shu (1977) and Whitworth and Summers 
(1985). As shown by these authors, there is a 2D con- 
tinuous set of solutions (taking into account the solutions 
which present weak discontinuities at the sonic point) de- 
termined for example by the value of the central density 
with bands of allowed and forbidden values. Two pecu- 
liar cases have been carefully studied, the so-called Larson- 
Penston and Shu solutions. The first case presents super- 
sonic velocities (up to 3.3 Cg for large radius, where Cg 
is the isothermal sound velocity) and is representative of 
very dynamical collapses. The second case assumes a qua- 
sistatic prestellar phase so that, att = 0, the density profile 
corresponds to the singular isothermal sphere (SIS) and is 
given by psis — cl/2TrGr^. A rarefaction wave which 
propagates outwards is launched and the collapse is inside- 
out. For both solutions the outer density profile is oc 
whereas in the neighbourhood of the central singularity, the 
density is cx r"-'^ '"'. 

Although the self-similar solutions depart significantly 
from the numerical calculations, they undoubtedly provide 
a physical hint on the collapse and the broad features de- 
scribed above appear to be generic and are observed in the 
simulations. Following the work of Foster and Chevalier 
(1993), various studies have focussed on the collapse of 
a nearly critical Bonnor-Ebert sphere (Ogino et ai, 1999; 
Hennebelle et ai, 2003). This scenario presents a number 
of interesting features which agree well with observations 
of dense cores like those observed in the Taurus molecular 



cloud (Tafalla et at., 1998; Bacmann et at., 2000; Belloche 
et at., 2002). Namely: (i) the density profile is approx- 
imately flat in the centre during the prestellar phase; (ii) 
during the prestellar phase there are (subsonic) inward ve- 
locities in the outer layers of the core, whilst the inner parts 
are still approximately at rest; (iii) there is an initial short 
phase of rapid accretion (notionaUy the Class phase), fol- 
lowed by a longer phase of slower accretion (the Class I 
phase). This last feature is an important difference with 
the self-similar solutions, which have a constant accretion 
rate. The typical accretion rates obtained numerically are 
between the value of the Shu solution (Ms is — c^/G) and 
the Larson-Penston solution (about 50 xcl/G). 

Motivated by the observations of much faster infall (see 
e.g. Di Francesco et ai, 2001), ttiggered collapses have 
been considered (Boss, 1995; Hennebelle etal, 2003, 2004; 
Motoyama and Yoshida, 2003). Much larger accretion rates, 
higher cloud densities and supersonic infall can be obtained 
in this context. A close comparison between a strongly trig- 
gered collapse model and the class-0 protostar 1RAS4A has 
been performed with success by Andre et al. (2004). 

3.1.2. Second Collapse and formation of a young stellar 
object 

When the density becomes larger than ~ 10^" cm~'^ the 
gas becomes optically thick. The isothermal phase ends and 
the thermal structure of the collapsing cloud is nearly adi- 
abatic. A thermally supported core forms {Larson, 1969; 
Masunaga et al., 1998). When matter piles up by accretion 
onto this hydrostatic core, its temperature and density in- 
crease because of the sttonger self-gravitating field. When 
the density of the first Larson core reaches about 10~^ g 
cm~"^, temperature is about 2000 K and the H2 molecules 
start to dissociate {Saumon et al, 1995). Most of the grav- 
itational energy goes into molecular dissociation energy so 
that the effective adiabatic exponent, 7 = 1+ ^"^"^ drops to 
about 1.1, significantly below the critical value 7=4/3 {Lar- 
son, 1969; Masunaga and Inutsuka, 2000). Thermal pres- 
sure is therefore unable to support the hydrostatic core and 
the collapse restarts. 

During the second collapse the temperature is roughly 
constant and close to 2000 K. When all the H2 molecules 
have been dissociated into atomic hydrogen, the effective 
adiabatic exponent increases again above 7=4/3 and the star 
forms. The timescale of the second collapse is about the 
freefall time of the first Larson core, ~ 1 yr, very small com- 
pared with the timescale of the first collapse which is about 
1 Myr 

Both the first and second Larson cores are bounded dur- 
ing all the collapse of the cloud by an accretion shock in 
which the kinetic energy of the infalling material is con- 
verted into heat. The effect of the accretion shock onto the 
protostar has been first considered by Stahler et al. (1980) 
and Stahler (1988). The influence of accretion on the evo- 
lution of the protostar will be examined in §4. 
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3.2. Influence of rotation and magnetic field 

Here we examine the main influence of rotation and 
magnetic field on the cloud collapse, leaving aside 3D ef- 
fects which are considered in §3.3. 

3.2.1. Effects of rotation 

Rotation induces a strong anisotropy in the cloud, slow- 
ing down and finally stopping the equatorial material. Ul- 
rich (1976) studied exact solutions for a rotating and col- 
lapsing cold gas and showed that the equatorial density of 
the collapsing envelope is larger than in the absence of ro- 
tation. This has been further confirmed by Terebey et al. 
(1984) using an analytical solution which generaUses the 
collapse of the SIS {Shu, 1977) in the case of a slowly ro- 
tating cloud. In the case of a 1 Mq initially slowly rotating 
core (/3 = Erot/Egrav - 2%), Hennebelle et al. (2004) 
estimate that the equatorial density of the collapsing enve- 
lope in the inner part of the cloud ( < 2000 AU) can be 2 to 
3 times higher than the axial one for a slow collapse and up 
to 10 times higher in case of strongly compressed clouds. 

The formation, growth and evolution of the rotationally 
supported disk has been modeled analytically by Cassen & 
Moosman (1981) and Stahler et al. (1994). The growth 
of the disk drastically depends on the angular momentum 
distribution, j. The centrifugal radius is about: ~ 
p /GMint where Mint is the mass inside the sphere of ra- 
dius rd- Therefore, for initial conditions corresponding to a 
SIS in soUd body rotation. Mint oc r and j oc r"^, implying 
ra oc Mf^^. On the contrary, starting with a uniform density 
sphere in solid body rotation. Mint oc and rd oc M^^^, 
which implies much bigger disks. Such disks are indeed 
found in hydrodynamical simulations of collapsing dense 
core initially in slow rotation. For 1 Mq dense cores with 
(3 ~ 2% the size of the disk during the class-0 phase is 
about 200 AU. 

The effect of the rotation on the forming protostar itself 
has been weakly explored. 2D equilibrium sequences of 
rotating protostars have been calculated by Durisen et al. 
(1989). 

3.2.2. Effects of magnetic field 

Magnetic field has been proposed to be the main support 
of the dense cores against the gravitational collapse {Shu 
et al., 1987) and the explanation for the low star forma- 
tion efficiency in the Galaxy. Although this theory is now 
challenged by the origin of the support being mainly due 
to turbulence (see Mac Low & Klessen, 2004 for a recent 
review), magnetic field certainly plays an important role in 
the formation of the protostar. 

The magnetically controlled collapse has been care- 
fully investigated with ID numerical simulations (e.g. 



Mouschovias et al, 1985). It has been found that the col- 
lapse proceeds in 2 main phases, first a quasi-static con- 
traction of the flattened cloud occurs through ambipolar 
diffusion and second, once a supercritical core has devel- 
oped, it collapses dynamically. Quantitative estimates of the 
prestellar cloud lifetime are given in Basu and Mouschovias 
(1995). In strongly subcritical clouds (initial mass-to-flux 
ratio over critical mass-to-flux ratio smaller than 1/10) the 
formation of the protostar requires about 15 freefall times 
whereas in a transcritical cloud (initial mass-to-flux ratio 
equal to critical mass-to-flux ratio), it requires about 3 
freefall times. Ciolek and Basu (2000) showed that the col- 
lapse of the well studied prestellar cloud, LI 544, is compat- 
ible with this core being transcritical. Note that, although 
the ambipolar diffusion time scale is much larger than the 
admitted star formation timescale, namely a few dynamical 
timescales, recent 2D simulations of compressible turbu- 
lence by Li and Nakamura (2004) suggest that enhanced 
ambipolar diffusion occurs through shock compression. 

The transfer of angular momentum is another important 
effect of magnetic fields. It occurs through the emission 
of torsional Alfven waves which carry away the angular 
momentum {Mouschovias and Paleologou, 1980; Basu and 
Mouschovias, 1995). Since this process is more efficient if 
the rotation axis is perpendicular to the field lines (instead 
of parallel), alignment between the magnetic field and the 
rotation axis is rapidly achieved. During the supercritical 
core formation epoch the angular velocity achieves a limit- 
ing profile proportional to l/r {Basu, 1998). Such a profile 
leads to centrifugal disks growing as oc Mint and thus 
intermediate between the very massive disks found in hy- 
drodynamical simulations and the low-mass disks predicted 
by the SIS in solid body rotation model. 

A very important difference between hydro and MHD 
cases is the presence of outflows in the latter ones, which 
have been found only recently in numerical simulations of 
collapsing protostellar core. They are described in the next 
section. 

Finally magnetic fields may induce a different mode of 
accretion. Motivated by the observations of T Tauri stars, 
which are surrounded by a disk from which they accrete ma- 
terial while having rotation velocities too small to be com- 
patible with the conservation of angular momentum, Konigl 
(1991) proposed that most of the accreted matter may be 
channeled along the magnetic field lines from the disk to the 
poles of the star. The angular momentum is then extracted 
from the infalling gas by the magnetic field. The accretion 
onto the star occurring over a small fraction of its surface, 
significant differences with the case of spherical accretion 
are expected {Hartmann et al., 1997), an issue addressed in 
§4. 

3.3. Tliree dimensional models 
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Fig. 2. — Accretion rate (in Mq/jt) and average angle of accretion during the 3D simulation of the second collapse of a 10 ^ Mq core. 
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Fig. 3. — Radial density and temperature profiles along the equatorial direction during the collapse of a 10~^ Mq core, for 4 different 
time steps. Left column: ID (spherical) collapse (dash-dot=2.860 yr, dash=2.866 yr, dot=2,879 yr, solid=2.906 yr); right colimm: 3D 
collapse of a rotating core (dash-dot=3.594 yr, dash=3.606 yr, dot=3.634 yr, solid=3.704 yr). Note the different behavior of the accretion 

shock in the two cases. 
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3.3.1. Axisymmetry breaking, transport of angular mo- 
mentum and fragmentation 

One of the main new effects which appear in 3D calcula- 
tions of a collapsing cloud is the axisymmetry breaking of 
the centrifugal disk. This occurs when its rotational energy 
reaches about 40% of its gravitational energy. Strong spiral 
modes develop which exert a gravitational torque leading to 
a very efficient outwards transport of angular momentum, 
allowing accretion onto the central object to continue. This 
effect has been modeled analytically (Laughlin & Rozyczka, 
1996) and found by many authors in the numerical simula- 
tions (e.g. Matsumoto and Hanawa, 2003). 

The fragmentation of the dense cores and the formation 
of multiple systems is one of the main challenges of the field 
and entire chapters of this book are dedicated to this subject. 
We refer to those as well as to the review of Bodenheimer 
et al. (2000b) for a comprehensive discussion of this topic. 

3.3.2. Multidimensional treatment of the second col- 
lapse 

The second collapse leading to the formation of the pro- 
tostar has been modeled in 2 or 3D by various authors with 
two main motivations, namely modelUng outflows and jets 
and explaining the formation of close binaries. Due to the 
large range of dynamical scales involved in the problem, the 
first calculations started from the first Larson core {Boss, 
1989; Bonnell and Bate, 1994). With the increase of com- 
putational power, calculations starting from prestellar core 
densities (e.g. 10* cm~^) have been performed {Bate, 1998; 
Tomisaka, 1998; Banerjee and Pudritz, 2005). For com- 
putational reasons, the radiative transfer has not been cal- 
culated self-consistently yet. Instead, piecewise polytropic 
equations of state which mimics the thermodynamics of the 
cloud are often used {Bate, 1998; Jappsen et al, 2005). 
More recently Banerjee and Pudritz (2005) used a tabulated 
cooling function which takes into account the microphysics 
of the gas with an approximated opacity. 

Bonnell and Bate (1994) conclude that fragmentation is 
possible during the second collapse. However since the 
mass of the stars is of the order of the Jeans mass, it is 
very small (0.01 Mq) and therefore they have to accrete 
most of their final mass. Banerjee and Pudritz (2005) form 
a close binary (with a separation of about 3 Rq) as well 
in their MHD adaptive mesh refinement calculations. Like 
Bate (1998) they find that inside the large outer disk (60- 
200 AU) an inner disk of about 1 AU forms. 

Tomisaka (1998) and Banerjee and Pudritz (2005) re- 
port outflows and jets during the collapse which contribute 
to carry away large amount of angular momentum. The 
physical mechanisms which is responsible for the launching 
of these outflows can be understood in terms of magnetic 
tower {Lynden-Bell, 2003). An annulus of highly wound 
magnetic field lines is created by the rotational motions and 



pushes the surrounding infalling material outwards. The 
physics involved in the jet is somehow different and based 
on the magneto centrifugal mechanism proposed by Bland- 
ford and Payne (1982). 

In the 3D simulations below, we investigate inner 
core formation resulting from the collapse of a 10~^ Mq 
Bonnor-Ebert sphere with densities and temperatures char- 
acteristic of the second core, namely p ~ 10~^ g.cm~^ and 
T ~ 2000 K. We focus on the influence of tri-dimensional 
effects on the accretion geometry and on the inner profile of 
the core. Fig. 2 displays the evolution of the accretion rate 
M during the second collapse as well as the average angle 
of accretion (cos 9), i.e. the average angle between the ver- 
tical axis in spherical coordinates and the infalling gas. As 
seen in the figure, the accretion rate decreases immediately 
from a large value close to the Larson-Penston prediction 
to a smaller Bondi-Hoyle or Shu like value, / G, and ac- 
cretion occurs over a very limited fraction of the protostar 
surface, {cos 6) < 0.3 (spherical accretion would imply 
(cos 9) = 0.5), so that most of the surface can radiate freely 
its energy. This is important for the subsequent evolution 
of the object, as examined in the next section. The con- 
sequences of 3D effects on the density and temperature 
profiles of the protostar are illustrated in Fig. 3 which dis- 
plays the equatorial density and temperature profiles of the 
second Larson core at 4 time steps. Rotation leads to lower 
central densities and temperatures and to a more extended 
central core, as noted already by Boss (1989). These fea- 
tures are relevant for the internal energy transport - radia- 
tion vs convection - and the initial deuterium burning. They 
also confirm that spherical collapse, although providing in- 
teresting qualitative information, cannot provide accurate 
initial conditions for PMS tracks as it will overestimate (i) 
the internal temperature of the protostar and (ii) the surface 
fraction covered by accretion, thus preventing the object to 
contract at a proper rate. 

4. EFFECT OF ACCRETION ON THE EARLY 
EVOLUTION OF LOW-MASS OBJECTS 



4.1. Observed accretion rates 

Intensive investigations of accretion in young clusters 
and star formation regions show signatures of this pro- 
cess over a wide range of masses, down to the substellar 
regime (see recent work by Kenyan et al, 2005; Mohanty 
et al., 2005; Muzerolle et al., 2005, and references therein). 
In the youngest observed star forming regions, such as p- 
Ophiuchus with an age < 1 Myr, the fraction of accre- 
tors is greater than 50%, independent of the mass {Mo- 
hanty et al, 2005). This fraction decreases significantly 
with age, a fact interpreted as a decrease of the accretion 
rates below the observational limits, < 10~^^M©/yr. The 
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timescale for accretion rates to drop below such a mea- 
surable limit is ^ 5 Myr. In some cases, however, ac- 
cretion continues up to ~ 10 Myr. Note, however, that 
these age estimates for young clusters remain very uncer- 
tain, since they are usually based on evolutionary tracks that 
are not rehable at such ages {Barajfe et al, 2002). Indeed, 
as demonstrated in Bamffe et al, (2002), unknown initial 
conditions and unknown convection efficiency (mimicked 
in stellar evolution calculations by the mixing length pa- 
rameters) during the early PMS contracting phase, charac- 
terized by short Kelvin-Helmholtz timescales ( < 10^ yr), 
can affect drastically the contraction track of a young ob- 
ject in the Herzsprung-Rusell (HR) diagram. Therefore the 
age and/or mass of young objects can not be determined 
accurately from observations, leading to very uncertain in- 
ferred disk lifetimes. However, even though the absolute 
timescales are uncertain, the trend of accretion rates de- 
creasing with time is less questionable. A sharp decrease 
of accretion rates with mass is also observed, with a cor- 
relation M cx M^, all the way from solar mass stars to 
the smallest observed accreting brown dwarfs, i.e ~ 0.015 
Mq (Muzerolle et al., 2005). Typically, in the low mass 
star regime (M ^ 0.2 — 1 Mq), the accretion rates vary 
between lO~^''M0/yr and lO^^M0/yr, whereas below ~ 
0.2 Mq and down to the BD regime, accretion rates range 
from 5 X 10"^ Mq/jt to l(r^^MQ/yr (Muzerolle et al, 
2003; Natta et al, 2004; Mohanty et al, 2005). Last but 
not least, observations now show similarities of accretion 
properties between higher mass stars and low mass objects, 
including brown dwarfs, suggesting that stars and brown 
dwarfs share similar formation histories. 

4.2. Modeling the effect of accretion in young objects 

On the theoretical front, Stahler (1983, 1988) has in- 
vestigated the effect of spherical accretion onto protostars, 
defining the concept of a birth line, a locus in the HR dia- 
gram where young stars first become optically visible when 
accretion ends. Stahler suggested that when the infall of 
material onto the protostar, responsible for its obscuration, 
ceases abruptly, the central object becomes an optically 
bright T Tauri star. 

Since this benchmark work, progress in the observations 
of young objects have now shown that accretion occurs 
rapidly through a disk, as discussed in §3.2 and 3.3 and il- 
lustrated in Fig. 2. The timescale for disk accretion is much 
longer than the strongly embedded protostellar phase, as il- 
lustrated by the short lifetime of the class-0 objects com- 
pared with class-1. Several studies have investigated the ef- 
fect of accretion geometry on evolutionary tracks for low- 
mass and high-mass stars. These calculations generally as- 
sume that (i) accretion takes place over a small fraction 5 
of the stellar surface and (ii) a dominant fraction of the ac- 
cretion luminosity is radiated away and thus does not mod- 
ify the protostar internal energy content {Mercer-Smith et 
al, 1984; Palla and Stahler, 1992; Hartman et al, 1997; 



Siess et al, 1997), in contrast to the assumptions of Stahler 
(1988). Under these conditions, the luminosity of the ac- 
creting object is given by: 

L = 6 ■ Lace + Ljj 

On the right hand side of Eq.5, the first term is the accreted 
luminosity, supposed to be entirely radiated away, Ld is 
the D-burning luminosity, including freshly accreted deu- 
terium, while the last term stems from the extra entropy at 
constant time due to the accreted mass (where m = m(m') 
is the accreted rate per mass shell). The first assumption (i) 
is indeed relevant for thin disk accretion from a boundary 
layer or for magnetospheric accretion where the gas falls 
onto the star following magnetic accretion columns. It im- 
plies that most of the stellar photosphere can radiate freely 
and is unaffected by a boundary layer or accretion shocks. 
The second assumption (ii) depends on the details of the ac- 
cretion process, which remain very uncertain. In a attempt 
to study the impact of such an assumption on evolutionary 
models, one can assume that some fraction of the accreted 
matter internal energy is transferred to the protostar outer 
layers, the other fraction being radiated away. This extra 
supply of internal energy, per unit mass of accreted matter, 
is proportional to the gravitational energy, eGM/R, with 
e < 1 a free parameter. As pointed out by Hartmann et al. 
(1997), the structure of an accreting object before or after 
ignition of deuterium, and the fact that it will be fully con- 
vective or will develop radiative layers, strongly depends 
on e, and to a lesser extent on assumption (i). For large 
values of e, convection can indeed be inhibited, even after 
deuterium ignition (see, e.g., Mercer-Smith et al., 1984). 
Deuterium burning in the protostellar phase is also a cen- 
tral issue. The key role played by deuterium burning on the 
properties of an accreting object and its location in the HR 
diagram was highhghted by Stahler (1983, 1988). Whether 
the deuterium fusion occurs in a fully convective object or 
in radiative layers is thus an important issue that affects sig- 
nificantly the structure of an accreting object. 

Assuming that only a very small fraction of the thermal 
energy released by accretion is added to the stellar interior, 
most of it being radiated away, Hartmann et al (1997) (see 
also Siess et al, 1999) showed that, depending on its evo- 
lutionary stage, an accreting low mass star expands less or 
contracts more than a non accreting similar object. Con- 
sequently, an accreting object looks older in a HR diagram, 
because of its smaller radiating surface for the same internal 
flux, compared to a non accreting object at the same mass 
and age. This stems essentially from the accretion timescale 
becoming of the order of the Kelvin-Helmholtz timescale, 
for a given accretion rate, M / M w Ikh, so that the con- 
tracting object does not have time to expand to the radius 
it would have in the absence of accretion. An extension of 
these studies to the brown dwarf regime confirms these re- 
sults, in the case of significant accretion rate and no thermal 
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energy addition due to accretion (e = 0) (Gallardo, Barajfe 
and Chabrier, in preparation). Fig. 4 shows the effect of ac- 
cretion on the radius of an object with initial mass 0.05 Mq, 
with accretion rate M = 10'^ Mq/jv and ^ = 0, e = 0. At 
any time, its structure is more compact than that of a non 
accreting object of same mass (dashed curve in Fig. 4), as 
mentioned above and as expected for accretion onto a fully 
convective object (Prialnik and Livio, 1985). The smaller 
radius, and thus the smaller luminosity, affect the location 
of the accreting brown dwarf in a HR diagram, as illustrated 
in Fig. 5. As seen on this figure, assigning an age or a mass 
to an observed young object of a given luminosity using 
non-accreUng tracks can significantly overestimate its age, 
at least with the present accretion parameters. The effect of 
various accretion rates (see below) and of finite values of e 
is under study. This again illustrates the uncertainty in age 
determination based on evolutionary tracks at young ages. 

4.3. Perspectives 

The calculations presented above for an accreting brown 
dwarf have been done with no or small transfer of internal 
energy from the accretion shock to the brown dwarf inte- 
rior. But our understanding of accretion mechanism is still 
too poor to exclude the release of a large amount of energy 
due to accretion at deep levels. As mentioned previously, 
although current observations indicate low accretion rates 
(M < 10^^ Mq/jt) for brown dwarfs at ages >1 Myr, 
they also point to rates decreasing with increasing time, 
suggesting significantly larger accretion rates at early times 
(^ 1 Myr). If large amounts of matter are accreted, even 
through a disk, one expects a significant amount of thermal 
energy to be added to the object internal energy (Hartman 
et al., 1997; Siess et al, 1999). In which case we expect 
important modifications of the structure of the surface lay- 
ers, with possible inhibition of convection as predicted for 
more massive objects {Mercer-Smith et al., 1984; Palla and 
Stahler, 1992), and thus a larger impact on ages and lo- 
cations in the HR diagram than displayed in Fig. 5. Such 
effects of accretion need to be explored in details in order to 
get a better characterization of their impact on the early evo- 
lution of low mass stars and brown dwarfs and thus of the 
uncertainties in mass and/or age determinations for young 
low mass objects. 

5. BROWN DWARF VS PLANET: OBSERVABLE 
SIGNATURES 

The "planetary status" of objects below the deuterium- 
burning Umit, ~ 13 Mjup {Saumon et al, 1996; Chabrier 
et al, 2000), remains the subject of heated debate. The 
debate was recently intensified by the direct image of an 
object below this mass limit, 2M1207b, orbiting a young 
brown dwarf at a projected orbital distance > 55 AU (Chau- 
vin et al, 2004). The present lAU working definition of a 



"planet" relies primarily on mass - not on the formation 
mechanism. However, to understand the formation mech- 
anisms of very low-mass objects, it is critical that we be 
able to single out those which formed in a disk by a three 
step process as described in §2.1 (core-accretion followed 
by gas-capture) from low-mass, no deuterium burning ob- 
jects which potentially formed by gravitational collapse of 
a molecular cloud fragment. According to the definition 
adopted in the present review, the former would be iden- 
tified as genuine planets while the latter would be brown 
dwarfs. It is interesting, by the way, to note that D-burning 
is advocated to distinguish BDs from planets, whereas stars 
with masses below and above the limit for ignition of the 
CNO cycle share the same "star" denomination. A com- 
mon "brown dwarf denomination should thus be used for 
D-buming or not D-buming BDs. Indeed, D-buming is 
essentially inconsequential for the long term evolution of 
these objects, in contrast to steady hydrogen burning which 
yields nuclear equilibrium and determines completely the 
fate of the object, star or brown dwarf (see e.g. Chabrier 
and Barajfe, 2000, Fig. 2 and 6). 

In the coming decades, direct imaging surveys are cer- 
tain to yield a sizeable number of objects below 13 Mjup 
orbiting stars and brown dwarfs beyond a few AU's - a 
region unlikely to be well sampled by radial velocity sur- 
veys. Without a disk signpost, it will be difficult to dis- 
tinguish long-period planets from very low mass brown 
dwarfs, based on their different formation history. A very 
low-mass brown dwarf (that never burned deuterium) could 
well be mistaken for a massive planet (see §2.2.1). Observ- 
able features that can distinguish between these two types 
of objects are greatly needed. 

Possible formation signatures could be contained in the 
atmospheric abundance patterns of planets and their mass- 
luminosity relationships. As mentioned in §2.1, a planet 
recently forged in a disk by the three-step process will ex- 
perience a brief period of bombardment which enriches its 
atmosphere and interior in metals compared to its parent 
star abundances, as observed for our jovian planets {Bar- 
shay and Lewis, 1978; Fegley and Ladders, 1994; Bezard 
et al, 2002, see also the chapter by Marley et al). Brown 
dwarfs, on the other hand, should retain the abundance pat- 
tern of the cloud from which they formed and, in the case 
of BDs in binaries, should have abundances similar to their 
primary star. The metalhcity distribution of planet-hosting 
stars found by radial velocity surveys already suggests that 
planet formation is favored in metal rich environments thus 
making an abundance test even more attractive. Recent in- 
terpretations of Spitzer observations for two extrasolar plan- 
ets, are suggestive of non-solar C and O abundances (see the 
chapter by Marley et al. and references therein). 

Enhanced metallicity leaves its mark on the interior, at- 
mospheric structure, and emergent spectrum in a variety 
of ways. As mentioned in §2.2.2, the presence of a large 
heavy element content in the planet interior will affect its 
mechanical structure, i.e. its mass-radius relationship. It 
will also modify its atmospheric structure. Fig. 6 compares 
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Fig. 6. — Left: Temperature versus pressure for a young Jupiter-mass planet atmosphere model with solar and 5 times solar metal 
abundances (i.e., [Fe/H]=0.7). Right: Model spectra for the same conditions. The Spitzer IRAC filter is indicated. 
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Fig. 4. — Effect of accretion on the mass-radius relationship of 
an accreting brown dwarf with initial mass 0.05 Mq and accretion 
rate M = lO~*M0/yr (solid line). The dashed line indicates the 
radius of a non-accreting object with same mass and same age as 
its accreting counterpart. Ages for the accreting object, in Myr, 
are indicated by the numbers. 



Fig. 5. — Evolution in the HR diagram of an accreting brown 
dwarf, with initial mass 0.05 Mq and accretion rate M = 
lO~**Af0/yr (solid line). The vertical dashed lines are cooling 
tracks of non-accreting low mass objects, with masses indicated 
near the curves (from 0.05 A/0 to 0.1 A/0). The square sym- 
bols indicate the position of non-accreting objects with the same 
age (indicated by the numbers near the squares, in Myr) and same 
mass as the accreting counterpart (indicated by a triangle just be- 
low the corresponding square). The numbers in brackets (close to 
the triangles) give the age (in Myr) of a non-accreting object at the 
position indicated by the triangle. 
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model atmospheric structures for a young (Tg// = 1200K, 
\og{g) = 4.0), cloud-free, non-irradiated planet mass ob- 
ject with solar and 5 times solar abundances. As the atmo- 
spheric opacities increase with increasing metalhcity, a nat- 
ural warming occurs in the deeper layers of the atmosphere. 
This warming of the atmospheric structure will have a di- 
rect impact on the evolution and predicted mass-luminosity 
relationship. 

The right panel of Fig. 6 illustrates the spectral differ- 
ences between these two models. Clearly the most promi- 
nent effect is seen around 4.5 /im where the increased ab- 
sorption is due to an increase in CO. Since this CO band 
falls in the Spitzer IRAC (3 to 8 /zm) coverage, significant 
metallicity enhancements in planets could set them apart 
from typical brown dwarfs on an IRAC color-color dia- 
gram. There is also a noticeable increase in the JsT-band 
(~ 2.2/xm) flux. 

The main purpose of this section is simply to point out 
one avenue to explore; however, clearly a great deal of work 
must be done before a concise picture of the expected abun- 
dance patterns in planets is developed. Non-equilibrium 
CO chemistry, for example, is predicted to occur in cool 
so-called T-dwarf BDs (Fegley and Ladders, 1996; Saumon 
et al., 2003). Moreover, brown dwarfs forming by gravi- 
tational collapse will certainly have abundance patterns as 
varied as their stellar associations, some being relatively 
metal rich, e.g., the Hyades {Taylor and Joner, 2005). Addi- 
tionally, metallicity effects in broad band photometry could 
well be obscured by other competing factors like gravity. 
Also, our own solar system planets show a range of C-O 
abundance ratios and varying levels of CO atmospheric en- 
hancement due to vertical mixing. Careful examinations 
of all these effects are necessary before any rehable spec- 
tral diagnostic can be used to distinguish low-mass brown 
dwarfs from planets. Such a diagnostic, however, has the 
virtue to rely on a physical distinction between two distinct 
populations in order to stop propagating confusion with im- 
properly used "planet" denominations. 

6. Conclusion 

In this review, we have explored (non exhaustively) our 
present understanding of the formation and the early evo- 
lution of gaseous planets and protostars and brown dwarfs. 
We now have consistent calculations between the planet for- 
mation, and thus its core mass and global heavy element 
enrichment, and the subsequent evolution after disk dissi- 
pation. These calculations are based on a revised version 
of the core accretion model for planet formation, which in- 
cludes planet migration and disk evolution, providing an ap- 
peahng scenario to solve the long standing timescale prob- 
lem in the standard core accretion scenario. Uncertainties 
in the initial conditions of planet formation, unfortunately, 
lead to large uncertainties in the initial radius of the new 
bom planet. Given the dependence of the thermal Kelvin- 



Helmholtz timescale on radius, this translates into large un- 
certainties on the characteristic luminosity of young plan- 
ets, over about 10^ yr for a 1 Mjup planet. Thus, it is im- 
possible to say whether young planets are bright or faint and 
what is their initial gravity for a given mass and therefore 
whether their evolution wiU differ from the one of young 
low-mass brown dwarfs. Conversely, future observations of 
young planets in disks of reasonably well determined ages 
will enable us to constrain these initial conditions. 

We have explored the effect of multidimensional col- 
lapse on the accretion properties and mechanical and ther- 
mal structures of protostellar cores. These calculations 
demonstrate that, within less than a free fall time, accre- 
tion occurs non-spherically, covering only a very limited 
fraction of the surface, so that most of the protostar sur- 
face can radiate freely into space. Spherical collapse is 
shown to overestimate the inner density and temperature of 
the presteUar core, yielding inaccurate initial conditions for 
PMS contracting tracks. This is important for initial deu- 
terium burning and for energy transport, 3D inner structures 
having cooler temperatures and more extended cores. This 
issue, however, can not be explored correctly with numeri- 
cal tools available today as it requires multidimensional im- 
plicit codes. The effect of accretion on the contraction of 
young brown dwarfs was also explored. Even though pre- 
liminary, these calculations confirm previous results forpre- 
main sequence stars, namely that, for accretion timescales 
comparable to the Kelvin-Helmholtz timescale, the accret- 
ing object has a smaller radius than its non accreting coun- 
terpart, for the same mass and age, and thus has a fainter 
luminosity. This smaller radius, along with the possible 
contribution from the accretion disk luminosity, can lead to 
inaccurate determinations of young object ages and masses 
from their location in an HR diagram, stressing further the 
questionable validity of mass-age calibrations and disk life- 
time estimates from effective temperature and luminosity 
determinations in young clusters. These calculations also 
suggest that, because of the highly non-spherical accretion, 
young stars or brown dwarfs will be visible shortly after the 
second collapse and, depending on their various accretion 
histories, will appear over an extended region of the HR 
diagram, even though being coeval. This seems to be sup- 
ported by the dispersion of low-mass objects observed in 
young stellar clusters or star forming regions when placed 
in an HR diagram (see e.g. Fig. 1 1 of Chabier, 2003). This 
suggests that, in spite of all its merits, the concept of a well 
defined birth line is not a correct representation, as star for- 
mation rather leads to a scatter over an extended area in the 
HR diagram. 

Finally, we suggest the deuterium-burning official dis- 
tinction between brown dwarfs and planets to be abandoned 
as it reUes on a stellar (in a generic sense, ie including 
brown dwarfs) quasistatic formation scenario which now 
seems to be superseded by a dynamical gravoturbulent pic- 
ture. Star formation and planet formation very Ukely over- 
lap in the ^ few M^^-p range and a physically motivated 
distinction between these two different populations should 
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reflect their different formation mechanisms. Within the 
general paradigm that brown dwarfs and stars form pre- 
dominantiy from the gravoturbulent collapse of a molecu- 
lar cloud and should retain the composition of the parent 
cloud and that planet form dominantly from planetesimal 
and gas accretion in a disk and thus should be significantly 
enriched in heavy elements compared to their parent star, 
we propose that these distinctions should be revealed by 
different mechanical (mass-radius) and spectroscopic sig- 
natures. Further exploration of this diagnostic is necessary 
and will hopefully be tested by direct obervations of gen- 
uine exoplanets. 
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